An observational test of the Vainshtein mechanism 
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Modified gravity theories capable of genuine self-acceleration typically invoke a galileon scalar 
which mediates a long range force, but is screened by the Vainshtein mechanism on small scales. In 
such theories, non-relativistic stars carry the full scalar charge (proportional to their mass), while 
black holes carry none. Thus, for a galaxy free-falling in some external gravitational field, its central 
massive black hole is expected to lag behind the stars. To look for this effect, and to distinguish 
it from other astrophysical effects, one can correlate the gravitational pull from the surrounding 
structure with the offset between the stellar center and the black hole. The expected offset depends 
on the central density of the galaxy, and ranges up to ~ 0.1 kpc for small galaxies. The observed 
offset in M87 cannot be explained by this effect unless the scalar force is significantly stronger than 
gravity. We also discuss the systematic offset of compact objects from the galactic plane as another 
possible signature. 
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There has been a lot of interest in theories of modi- 
fied gravity that might explain the observed accelerated 
expansion of the universe [l|, • Theories capable of gen- 
uine self-acceleration - without invoking vacuum energy 
in the Einstein frame - are especially interesting, and 
they generally involve introducing a scalar (ip) that re- 
spects the galileon symmetry p — > ip + b + c^x^, where b 
is a constant and c p is a constant vector 4]. This scalar 
mediates a long-ranged force, the so-called fifth force in 
addition to the usual gravitational force between objects. 
Thanks to the Vainshtein mechanism Q, the scalar is 
screened on small scales, so that solar system constraints 
are satisfied. To see how it works, let us illustrate with 
the simplest galileon model, inspired by DGP 0, 0] ; the 
equation for the scalar (in Einstein frame) is 1 2 : 

U ^+^ [( D( P) 2 ~ dudutpWartp] = -SrraGT/ , (1) 

where is the trace of the matter energy-momentum, 
which for our purpose can be equated with (negative of) 
the matter density p, assuming it is non-relativistic. The 
constant a quantifies the scalar-matter coupling, and is 
generically of order unity, i.e. of gravitational strength. 
For instance, it takes the value l/y/6 in massive gravity 
models [§| . The mass scale m is generally of the order of 
the Hubble constant today m ~ Hq. We are interested in 
solutions of this equation in the quasi-static limit, mean- 
ing time derivatives can be ignored and □ — > V 2 . On 
large scales, the linear term V 2 ip dominates; Eq. (JTJ re- 
sembles the Poisson equation, with <p playing the role of 



the gravitational potential. A localized source p yields 
a profile tp that scales inversely with distance r. How- 
ever, as one approaches the source, the interaction term 
(second term on the l.h.s.) dominates, and simple power 
counting reveals p oc *Jr. Thus, at small distances, p 
is screened relative to the normal 1/r gravitational po- 
tential. The transition scale is known as the Vainshtein 
radius, and is roughly given by (GM/m 2 ) 1 ^ 3 for spher- 
ically symmetric configurations, where M is the mass 
of the source. For instance, the entire solar system fits 
within the Vainshtein radius of the Sun, about 0.1 kpc, 
greatly suppressing the scalar force sourced by the Sun. 
What makes the galileon model attractive from this view- 
point, is that it is the same nonlinear interaction that is 
responsible both for Vainshtein screening, and for self- 
acceleration 

An important property of Eq. (fT]) is that it can be 



rewritten in the form 9 M J M = 



linear function of derivatives of <p 



-T M M , where J M is a non- 



define a scalar "charge" Q = — J d 3 xT^ 



One can thus 
/ d 3 xp, 3 

which is none other than the mass M. 4 As shown 
in [{|, the scalar charge also quantifies the response of 
an object to an external field, i.e. an external gravi- 
tational $ cxt + scalar p ey ± field exerts a net force of 
Mx = — MV$ e xt — aQVp cyL t on an object of mass M 
and charge Q. It is worth emphasizing that the Vain- 
shtein mechanism does not suppress the scalar charge Q 
at all - indeed, the fact that Q = M enforces the equiv- 
alence principle, i.e. making the motion of an object in- 
dependent of its mass. What the Vainshtein mechanism 
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1 Here, ip is related to the standard notation it for the galileon by 
7r = cap. 

2 Two additional galileon symmetric interaction terms can be writ- 
ten down in the equation of motion: (Dip) 3 — 3D(/p(9 fl 9 I /(/p) 2 + 
2{d fi d^) 3 and {Dip) 4 - 6(n V ) 2 {d^d^f + 8D(p(d^ip) 3 + 
3[(c9 M <9„i/>) 2 ] 2 — 6(9 A i9^(/p) 4 . All results in this paper apply in 
the presence of any combination of the galileon terms. 



3 One might wonder if — J dPxT^ = j (fix p is too strong an as- 
sumption, since an apparently non-relativistic object might have 
relativistic components, e.g. gluons in protons. It turns out that 
as long as the object in question is stationary, and has a finite 
extent, - / d^xT^ = - f d 3 xT ° = f d 3 xp holds, by virtue of 
a tensorial virial theorem Hot . For further discussions, including 
the renormalization of Q by quantum effects, see [I0t lll|l . 

4 That is, aside from an exception described below. 
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suppresses is the external scalar field ip cx t that an object 
senses when it is close to the source of that field. 

The one exception to Q = M is compact objects, such 
as black holes. These are objects whose mass receives 
a significant contribution from the gravitational binding 
energy, a contribution that is not included in the scalar 
charge Q. Compact objects such as neutron stars thus 
have Q/M < 1, with black holes as the extreme limit 
where Q/M — 0. This is consistent with the notion of 
black holes having no hair, more specifically no galileon 
hair which we prove in a separate paper [12J (see also 
[IH EH)- 5 Thus, under some external <E> ext and </? C xt 
fields, a normal non-compact star (as well as dark matter 
particles) would fall according to x = — V<I>ext — aV</? C xt, 
while a black hole would fall according to x = — V$ cx t, 
insensitive to the scalar force. 6 

The challenge is to find situations where the scalar (p cxt 
is not already suppressed by the Vainshtein mechanism. 
Since both black holes and stars typically reside within 
galaxies whose Vainshtein radii greatly exceed their sizes, 
it would seem hopeless to observe the purported differ- 
ence in the rate of fall between a normal star and a black 
hole. 

The galileon symmetry helps save the day. The sym- 
metry tells us that given any solution to the scalar Eq. 
(TT]), one can always add a linear gradient, i.e. tp cxt with a 
constant V(^ ox t, and obtain another solution. 7 For any 
given object, whether such a linear gradient is present 
and how large it is, depends on the boundary conditions. 
For a galaxy, the boundary conditions are supplied by the 
surrounding large scale structure. Interestingly, as is re- 
cently demonstrated in a series of numerical simulations 
[17H2G ] , the galileon scalar indeed obeys linear dynamics 
on sufficiently large scales (> 10 Mpc), meaning the large 
scale structure generates a galileon field that is unsup- 
pressed by the Vainshtein mechanism 8 . The large-scale- 



If one thinks of black holes as vacuum solutions, having no scalar 
charge is certainly a solution, but the point of no-hair theorem is 
to show that it is in fact the only solution, and thus the collapse 
of an actual star would presumably lead to such a configuration. 
We should stress that our proof in [lT3| concerns only spherically 
symmetric black holes. We take it as suggestive that rotating 
black holes likely share the same no-galileon-hair property, but it 
remains to be proven. In any case, the argument for a suppressed 
Q/M for compact objects is quite robust. 

The fact that compact objects have a suppressed scalar charge 
is not new 15, ljj. What is new with the advent of recent mod- 
ified gravity models is the screening mechanism, which revives 
scalar-tensor theories that are otherwise already ruled out by 
solar system tests. 

This argument relies heavily on the galileon symmetry. It is an 
interesting open question whether the offset effect between stars 
and black holes would be observable in a theory like P(X), which 
exhibits Vainshtein screening without the galileon symmetry. 
In other words, while the Vainshtein mechanism does operate on 
small scales, the Vainshtein zones of nonlinear objects do not per- 
colate the universe. Note the Vainshtein radius, estimated from 
an isolated spherically symmetric object, can be a misleading 



structure-generated scalar field has a long wavelength, 
and can be approximated as a linear gradient on the 
scale of a galaxy. This linear gradient penetrates the 
Vainshtein zone of the galaxy, and can act unsuppressed 
on the galaxy and its constituents. In other words, the 
galaxy falls according to this unsuppressed scalar ip cxt in- 
duced by large scale structure. So do its constituent dark 
matter, non-compact stars, but not its compact objects. 
The most readily observable compact object is the cen- 
tral massive black hole if there is nuclear activity. The 
central black hole, lacking a scalar charge, does not re- 
spond to the scalar force, while the stars (and dark mat- 
ter particles) do. The net effect is that the black hole will 
lag behind the stars in their overall large-scale-structure- 
induced motion. In other words, the black hole will be 
offset from the center of the galaxy, or more precisely, 
from the minimum of the galactic gravitational poten- 
tial. The non-zero offset means there's an extra (purely 
gravitational, not scalar) tug on the black hole from the 
central region of the galaxy. This suffices to compen- 
sate for the lack of a scalar force on the black hole, and 
keep the black hole and stars in equilibrium, moving in 
tandem. One can estimate the size of the offset r by 
equating the extra scalar acceleration sensed by the stars 
a|V(/? 0X t| with the extra gravitational tug on the black 
hole GAf ga i(< r)/r 2 , where Af ga i(< r) is the mass en- 
closed within radius r of the galaxy. We find a displace- 
ment 

r = 0.1 kpc ( 2 -^\ { lg*Sgl ) f°-° 1M ° P H .(2) 

Here, we estimate ^ext by 2a$ ext , since the linear scalar 
ipext satisfies the same Poisson equation for the gravita- 
tional potential $ cxt , but with the source term scaled up 
by 2a (see Eq. [JJ). The typical gravitational accelera- 
tion |V$ cxt | can be estimated by the typical peculiar mo- 
tion multiplied by Hubble: 300km/s x 70km/s/Mpc ~ 
20 ( km/s) 2 /kpc. A more careful calculation of the rms 
|V$cxt| using the observed matter power spectrum gives 
a number fairly close to this (e.g. [21]). It should be 
kept in mind however that |V$ e xt| is a stochastic quan- 
tity, and its value depends on environment. The central 
density of po ~ 0.01M©/pc 3 is appropriate for dwarf or 
low surface brightness galaxies, where the effect is the 
largest. We provide explicit scaling with the relevant pa- 
rameters in Eq. ([2]) so that one can easily extrapolate to 
other values. 

In modeling the central region of a galaxy, several is- 
sues should be kept in mind. First of all, the above es- 
timate assumes the density profile is approximately flat 
close to the center. For fixed density at the core radius, 



concept when applied to more complex situations. The tensor 
structure of Eq. JTJ is such that screening works very differently 
in non-spherically symmetric situations. 
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a steep profile would imply a small offset. Fortunately, 
low central density galaxies where the offset is the largest 
also have fairly flat profiles e.g. 22]. Second, the rele- 
vant central density is the one averaged outside the black 
hole's sphere of influence. Materials within the sphere of 
influence would simply move with the black hole. It is the 
materials outside that are important for determining the 
offset. Thus, Bahcall-Wolf type cusps are not relevant 
for our considerations [Hj|. 

Galaxies for which both the stars and the central mas- 
sive black hole are readily observable are those with a low 
level nuclear activity, namely Seyferts. Our estimate in 
Eq. ([2]) suggests that the offset would be observable pref- 
erentially in small galaxies. Depending on the distance 
(up to say tens of Mpc), galaxies with a central den- 
sity in the range ~ 0.003 — 0.03 M Q /pc 3 , corresponding 
roughly to rotational velocities around ~ 30 — 120 km/s 
[22}, should have a measurable offset. Until recently, such 
small Seyferts have not been well explored observation- 
ally. A classic case of a dwarf galaxy containing an active 
nucleus is NGC 4395 [H. A good size sample (~ 30) of 
small Seyferts was recently reported by |25|. 9 

An interesting question is how fast the black hole 
moves relative to the stars, on its way to the equilib- 
rium (offset) position. Using the same set of parameters 
displayed in Eq. , we find a fairly small velocity of ~ 2 
km/s. At such a velocity, dynamical friction is unimpor- 
tant. The time it takes for the black hole to traverse 
the requisite distance is about 10 s years. This is fairly 
close to some estimates of the nuclear activity time-scale 
[26l ] , though the conditions for nuclear activity are rather 
uncertain. 

M87 is an interesting case, where its massive black hole 
is known to be offset from the (bulge) stellar center by 
a projected distance of ~ 7 pc 27]. Its central density 
is about po ~ 20M Q /pc 3 . For our effect to reproduce 
such an offset thus requires a ~ 8, assuming the external 
gravitational acceleration on M87 is typical. This is a 
scalar-matter coupling that is quite a bit stronger than 
gravitational. The more plausible explanations for the 
observed offset are: acceleration by an asymmetric jet, 
and gravitational wave recoil from a merger |27j |. Two 
other explanations considered are: the active black hole 
being a member of a binary, and Brownian motion. The 
former can be constrained by the lack of a jet precession, 
and can also be tested by monitoring the system over 
time. The latter is a small effect, giving an offset < 0.1 
pc in the case of M87 [27) . 

These alternative mechanisms raise a practical ques- 
tion: in the event one observes a black hole offset in 
a lower central density galaxy, consistent with a ~ 1, 
how does one disentangle the modified gravity effect from 
other astrophysical effects? One can exploit a distin- 
guishing feature of the modified gravity effect, that is, 
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the offset should be correlated with the gravitational pull 
of the surrounding large scale structure, in both its di- 
rection and strength. For instance, galaxies are expected 
to stream out of voids, their resident black holes should 
lag behind the stars in that streaming motion. Voids are 
especially interesting places to look because the scalar 
field is expected to be unscreened there 10 . However, we 
expect large scale structure to source a linear (i.e. unsup- 
pressed) scalar even away from voids. It would be very 
useful to map out the precise gravitational and scalar 
fields in our neighborhood, using fairly detailed knowl- 
edge of the mass distribution of the local universe |28| . 
The two large-scale-structure-generated fields <J> ex t and 
ifcxt are expected to roughly align, but an accurate map 
would aid in isolating the modified gravity effect from 
other astrophysical effects. Note also that typical astro- 
physical effects produce a velocity offset (~ 10 — 1000 
km/s) that is quite a bit larger than what modified grav- 
ity predicts. 11 

It is worth emphasizing that this offset effect - namely 
the lagging of compact objects behind stars in the over- 
all streaming motion of the host galaxy - is not confined 
to the central massive black hole. Any compact objects, 
regardless of its mass, will display this effect, though the 
effect is larger for a black hole than a neutron star i.e. 
the offset is expected to scale as 2GM/R, where R is 
the radius of the object. The advantage of the central 
massive black hole is that it is readily observable even 
in a distant galaxy, provided it is active. For other less 
massive compact objects, the best bet is to look within 
our own galaxy. One possible signature is to see if com- 
pact objects are systematically offset from the galactic 
plane (defined by the stars), in the opposite direction 
of our galaxy's streaming motion. Using numbers from 
Eq. ([2]), and adopting the solar neighborhood value of 
Po ~ 0.18 M Q / pc 3 , we estimate the offset to be about 2 
pc for black holes. It would be useful to compute this off- 
set more carefully by calculating the precise large-scale- 
structure-generated scalar force on the Milky Way. 

A natural question is: should we expect the same off- 
set effect for other screening mechanisms? Theories that 



A rough argument goes as follows. Eq. JTJ can be rewritten 
schematically in the form: Hg 2 d 2 ip+(HQ 2 d 2 ip) 2 ~ctp/p. Thus, 
voids where p/p < 1/a are natural places where one can consis- 
tently ignore the nonlinear term compared to the linear term on 
the left. However, we expect the linear/unsuppressed scalar to 
be present even away from voids, i.e. large scale structure gen- 
erates a scalar field with long wavelengths that can penetrate 
Vainshtein zones around nonlinear objects. 

An exception is Brownian motion which produces small 
velocity perturbations to the massive black hole: ~ 
UstarCMstar/Mj, h J 1 / 2 . Using the relation mass of black hole 
M h . h . ~ 10 812 M Q (u s tar/200km/s) 4 - 24 [H, the smallest galax- 
ies («star ~ 30 km/s) have the largest velocity perturbation 
of ~ 0.2 km/s. For such galaxies, the spatial offset is ~ 
»*core(A4star/M] 3 .h.) 1 ' /2 ~ 0.03 kpc [27ll , thus comparable to the 
modified gravity effect. Hence, it is important to use the corre- 
lation with large scale structure as a discriminant. 
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screen by means of scalar self-interactions of the poten- 
tial type, such as the chameleon [30] or the symmetron 
[3lj | . operate very differently from those that screen by 
derivative interactions i.e. Vainshtein. Non-Vainshtcin 
theories do not have a conserved scalar charge Q. In- 
deed, the scalar force on/from the Sun is screened in 
these theories by virtue of the Sun's suppressed Q/M. 
In chameleon and symmetron theories, objects with a 
gravitational potential similar to the Sun's (~ 10~ 6 ), 
or deeper, have Q/M ~ 0. Thus, main sequence stars 
and black holes fall at the same rate, and one expects 
no offset between the two. However, dark matter and 
stars can fall differently if the host galaxy has a suffi- 
ciently shallow gravitational potential This can lead 
to an interesting warping of the galactic disk, pointed 
out by Jain and VanderPlas [32[ . For other observational 
tests of chameleon/symmetron screening, see for instance 

H MM MM ■ 

In summary, we propose a test of Vainshtein screen- 



ing in galileon theories by comparing the rate of fall for 
compact objects versus non-relativistic stars. A positive 
detection of a difference will be a great boost to some 
of the recent ideas of modifying gravity to explain cos- 
mic acceleration. A negative detection can be used to 
put an upper limit on the scalar-matter coupling a. A 
bound reaching 0.1 is conceivable with existing data, and 
will suffice to rule out most of these recent ideas. Per- 
haps the most interesting observation is that relatively 
small scale, local data can shed light on the dark energy 
problem and the nature of gravity. 
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